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ABSTRACT 

The S Stars in the Galactic-center region are found to be on near-perfect Keplerian orbits around 
presumably a supermassive black hole, with periods of 15-50 yr. Since these stars reach a few percent 
of light speed at pericenter, various relativistic effects are expected, and have been discussed in the 
literature. We argue that an elegant test of the Einstein equivalence principle should be possible 
with existing instruments, through spectroscopic monitoring of an S star concentrated during the 
months around pericenter, supplemented with an already-adequate astrometric determination of the 
inclination. In essence, the spectrum of an S star can be considered a heterogeneous ensemble of 
clocks in a freely-falling frame, which near pericenter is moving at relativistic speeds. 
Subject headings: Galaxy: center, Relativistic processes 



1. INTRODUCTION 

The equivalence of gravitational and inertial mass, or 
that gravity can be cancelled by transforming to a freely- 
falling frame, was tested within classical celestial me- 
chanics to high precision by the end of the nineteenth 
century. After all, if Mercury had a gravitational con- 
stant differing from (say) Jupiter's by one part per mil- 
lion. Mercury's orbit would not have fitted classical dy- 
namics well enough to highlight the anomalous precession 
(43" per century, or 10~^) that was later explained by 
general relativity. 

The Einstein equivalence principle (EEP) adds to the 
classical or weak equivalence principle the further phys- 
ical postulate that special relativity holds locally in a 
freely-falling frame. A consequence of the EEP is that 
the effects of gravity on test particles are fully described 
by endowing spacetime with a Riemannian metric, and 
having the particles follow geodesies of that metric. A 
further consequence of the EEP is that the temporal met- 
ric component is, to leading order, given by gtt — 1-I-2"I>, 
I where 4> is the Newtonian potential. Hence gravity 
causes time dilation. 

Different aspects of the EEP have been verifie d by 
I multip le experiments, as reviewed extensively by iWilll 
(|2006D . In particular, null-redshift experiments test that 
gravity — whatever it does — does the same to clocks of 
different types. Pound-Rebka experiments verify that, 
in a static gravitational field, time is dilated by a fac- 

— 1/2 

tor ~ 1 — $. GPS satellites and binary pulsars 

are effectively freely-falling clocks moving alQ ^ 10^^ 
and ^ 10^^ respectively, and exhibit time dilation of 

9u' - 1 - 2$. 

But excellent as the existing experiments on the EEP 
are, it would be even nicer to have a laboratory with ex- 
periments on multiple materials, the whole freely-falling 
at relativistic speeds. In this Letter, we suggest that the 
S stars in the Galactic-center region approximate such a 
laboratory. Stellar spectra contain absorption lines from 
different atomic species, which can be regarded as differ- 
ent clocks. The EEP asserts that special relativity holds 

^ We use geometrized units: GM = c = 1, where M is the 
black-hole mass. 



in the star's frame for all atomic processes (Local Lorentz 
Invariance), irrespective of where the star is (Local Posi- 
tion Invariance), and these can be tested by observations 
of a single star. The remaining ingredient of the EEP, 
namely that a freely falling frame does not itself depend 
on the composition of the star (Universality of Free Fall) 
requires multiple stars. The nature of a violation of the 
universality of free fall is difficult to speculate on, how- 
ever one could imagine for example, fits from different 
S Stars' redshift curves all yielding different values for 
the black hole mass. Such a deduction would signal a 
violation of the UFF. 

Going beyond the EEP, yet another possibility is that 
a star with a relativistically significant gravitational self 
energy, i.e., a pulsar, may give different results. Such 
a result would be a violation of the Strong Equiva- 
lence Principle. It is expected that a significant num- 
ber of pulsars inhabit the inner milliparsec, and al- 
though none h ave yet been foun d, searches are currently 
iinderwav(Ma cquart et al]l201Cll ). 

The S stars achieve the highest speeds of any known 
geodesies. For example, S2 reaches v > 0.03 at pericen- 
ter, which it last passed in 2002 and will again in 2018. 
The pericenter speed is an order of magnitude faster than 
for known binary pulsars. The spatial scales at pericenter 
are of order a light day. Hence, even if the quantitative 
constraints are initially modest, the S stars would test 
the EEP at velocity and spatial scales not reached by 
other experiments. 

Astrometric and redshift observations of the S stars 
show them to be on orbi ts (so far) indistinguishable from 
pure Kepleri an (Gilless en et al.l 120091 : IChez et al.l 120081 : 
Gilles sen et al.l I2009al) . Keplerian elements are known 
to high accuracy for many stars. For example, for S2 
the orbital inclination is measured to be 135 ± 1°. Ex- 
pected relativistic effects have been discussed extensively 
in the literature. T hese include not only time dilation 
(Zuckcr et al. 2006), but also secular orbit precession 
(Rubilar & Eckart 2001; Wil l20Q8,) plus Newtoni an per- 
turbations from other stars (jMerritt et al.ll2010D . kine- 
matic effects due to space curvature and fram e-dragging 
(jKannan fc Sahal 120091: iPreto fc Sahal |2009[) . and com- 
posite redshift perturbations including light-path effects 
(jAngelil fc Salia.2010: Angelil et al. 2010.) . 
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The EEP implies a time dilation of 1 + 2/r. Hence a 
spectral line intrinsically at will be redshifted to v, 
where 



where 



= (1 



Via 



1 



(1) 



Redshift is conventionally defined as j/q/i^ — 1. It is how- 
ever, somewhat more convenient if one defines ln(i/o/z^) 
as the redshift, and we will do so. In any case, the 2/r 
term constitutes a redshift perturbation of O (w^) and is 
the strongest relativistic perturbation. In other words, 
the leading-order perturbation due to relativity is time 
dilation on a Keplerian orbit. Space curvature perturbs 
the redshift at O (u'^) , frame dragging due to black-hole 
spin at O (v^) , with additional physical effects continu- 
ing at higher orders (jAngelil et al.ll20l"ol ). One way to 
test the EEP would be to replace 2/r in Equation ^ by 
a/r and fit observed redshift curves for the para meter a. 
In fa ct, this has already been attempted (Zuck e~et al.l 
l2006f ) but the available data appeared not yet sufficiently 
accurate to overcome systematic uncertainties. We argue 
below that observations concentrated near pericenter will 
be particularly useful. 

Relativistic perturbations to the S stars increase with 
decreasing pericenter distance. In orbit fitting from spec- 
troscopy, as one looks to stars with longer periods, the 
perturbations to the redshift from relativity become over- 
whelmed by perturbations due to massive perturbers. Of 
all the compact objects orbiting the central black hole, we 
see only the brightest stars. Thousands of solar masses 
worth of compact objects and a possibly significant dark 
matter component is expected to inhabit the inner arc- 
seconds. Our lack of understanding of the form of this 
extended system may play a significant role in obscur- 
ing relativistic perturbations for large orbits. However, 
current models for the Newtonian perturbations are esti- 
mated to most likely be weaker than the time dilat ion ef- 
fect for the shortest-perio d stars such as S2(Schode l et al.l 
120091: iMerritt et all (20101) . Hence we wih disregard the 
the extended-mass component in our analysis. 

2. OBSERVABLES 

Consider a star on a pure Keplerian orbit. Elementary 
celestial mechanics gives the position of the star paramet- 
rically, that is, both coordinates and time are expressed 
as functions of the so-called eccentric anomaly tp. At 
time 

t= ^(^-esini/;) (2) 

the position of the star in its orbital plane is 

(x,?/, 2) = 1= \ cos -0 — e, \/l — sin -0, ) . (3) 

Here P is the period, e is the eccentricity, while a is the 
semimajor axis in units of the gravitational radius of the 
black hole. 

We now rotate the coordinate system, first by the ar- 
gument of periapsis w about the z axis, then by the in- 
clination / about the new x axis. The observer is now 
along the new z axis. The redshift is 



vr\ dz a 
In ^ = — + - 
1/ at r 



Ac f{e,iv,'tlj)+Aii 
1 — e cos V' 



(4) 



/(e, uj, tp) — \J \ — cosoj cos?/' — sinw sin?/' 
and we have introduced the coefficients 



sin/ a 
Ac = —;=^ and Ar = — . 

J a a 



(5) 
(6) 



which can be interpreted as the amplitudes of the clas- 
sical and relativistic contributions to the redshift. We 
also note that t above is the coordinate time of emission. 
The time of observation is of course t plus the light travel 
time, whose varying part is the R0mer time delay, which 
is the z coordinate of the star 



Ar 



p 
'2^ 



(cos 'ip — e) sin u + \/\ — sin cos uj 



(7) 



and can also be interpreted as the time-integral of the 
redshift. 

Classically, the redshift curve determines only the 
combination sin//-y/a, leaving the inclination unknown. 
Hence, for a given redshift curve, the inferred orbital 
speed becomes infinite if the orbit is face-on. Not sur- 
prisingly, relativity prevents that happening. From Q 
and (El) it follows that 



^2 

sin^ / = -r^a. 
Ar 



(8) 



Hence time-dilation breaks the inclination degeneracy if 
a is known, or allows a to be measured if / is known 
from astrometry. 

It is worth mentioning that another way in which rel- 
ativity can break the inclination degeneracy is through 
the well-known precession of 



Aw 



a(l-e2) 



(9) 



per orbit. For binary pulsars, the cumulative w pre- 
cession amounts to several degrees per year, allowing a 
and hence / to be inferred ( Brumberg et al.l 11975:) . For 
Galactic-center stars, however, the a; precession builds 
up much more slowly. Over a single orbit the redshift 
due to precession is '-^ O (u'^) because Ao; ^ O (w^) over 
one orbit. 

3. PARAMETER RECOVERY 
We treat the following seven parameters as unknown. 

1. The period P . 

2. An additive constant on t. 

3. The eccentricity e. 

4. The argument of pericenter w, which along with the 
inclination / sets the orbit orientation with respect 
to the observer. Since the observer is on the z axis, 
the nodal angle f2 corresponds to a rotation around 
the line of sight, which leaves the redshift curve 
invariant. 

5. The intrinsic frequency vq. This is the absolute line 
calibration, plus any offset due to the observer's ra- 
dial velocity. Changing z/q shifts the redshift curve 
vertically. The tangential motion of the observer 
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is neglected here, since the Sun's orbital speed of 
~ 200 km s^^ in the Galaxy would contribute no 
significant Doppler shift. 

6. A classical redshift amplitude Ac- 

7. A relativistic redshift amplitude Af;. In fact, a, I, a 
form a degenerate trio and are absorbed into the 
non-degenerate dummy parameters Ac and Aji via 



Astrometry independently measures the first four of 
these, as well as /. Astrometric observations also involve 
six other parameters: the position and proper motion of 
the Galactic center on the sky, the distance to the Galac- 
tic center, and fi. The mass of the black hole is not an 
independent parameter, since it is a function of P and a. 
For our purposes, only / from astrometry is essential. In 
testing the EEP, while we shall argue that observation of 
the star over a short time around pericenter suffices for 
spectroscopy, recovery of the inclination via astrometry 
can be done anywhere on the orbit. 

It is important to note vq and v need not refer to 
a single spectral lin e. The spectra of S stars (see e.g., 
iMartins et al.l l2008') contain multiple features. Most are 
early-type stars with H and He features, while about 10% 
are late- type stars with molecular and metal bands/lines 
and little or no H or He. Hence, S-star spectra could, 
in principle, test the equivalence principle for multiple 
atomic processes. If the stellar atmosphere does not 
change appreciably over an orbit, an observed spectrum 
can be cross-correlated on a logarithmic wavelength scale 
with a spectrum observed at some other epoch, and 
the cross-correlation peak would directly give the red- 
shift ln(z^o/'^) with vq an unknown constant. If different 
atomic/molecular species behave differently in a freely- 
falling frame, the shape of the cross-correlation curve 
would change. Alternatively, multiple spectral features 
could be fitted simultaneously with variable redshift. We 
do not, however, attempt to model the observable spec- 
tra explicitly in this paper. 

We now simulate the recovery of a as follows. We gen- 
erate 10 mock redshift data points of S2 taken over two 
months at pericenter, plus four additional data points, at 
±1, ±2yr around pericenter. The data are generate d with 
a = 2, and orbital parameters taken from Gillessc n et al.l 
(12009a). To them we add gaussian random noise at a dis- 
persion of 10 km s~^, and then fit via the seven parame- 
ters. We then assume / has been measured by astrometry 
and use ^ to recover a. Figure ([l} shows an example 
for a few mock data realizations at a fixed accuracy, and 
figure shows the dependency of the recovered value 
of alpha with the data accuracy. 

4. DISCUSSION 

Testing the equivalence principle using a combination 
of spectroscopy and astrometry seems possible in the 
near future. In comparing the spread in A^/An from 
mock data to the recovered value for / from real astro- 
metric data (illustrated in Fig. [ij, in testing the Equiv- 
alence Principle using S2, the current accuracy available 
from astrometry sits at a comfortable level. Spectro- 
scopic accuracy of S2 at 10 km s~^ is not yet available, 
but seems plausible with future observations. For the 
late- type star S35, which has a more favorable spectrum. 
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Fig. 1. — Recovery of a for ten mock data sets. The ratio A^/A^ 
is the slope of the hnes in the above plot, and is recovered from 
mock spectroscopic data of 14 data points with accuracy 10 km s~^ , 
concentrated around pericenter. The horizontal lines are the up- 
per and lower confidenc e levels for the recover ed inclination from 
astrometry, taken from [HiUessen et al.l l l2009al V The intersection 
point corresponds to the value of a for which both datatypes agree 
on the inclination. 
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Fig. 2. — The Itr level in the recovered value of a for spectroscopic 
data of different accuracy. At each data accuracy level, we have 
performed the same procedure illustrated in Figure [T] except that 
60,000 mock data realizations have been used. 

a fit error of 10 km s^^ has been achieved (|Gillessen et al.l 
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l2009a[) . We remark that any systematic errors that do 
not change between observations are harmlessly absorbed 
into vq. 

Naturally, in detecting relativistic effects on S Stars, 
data during pericenter passage are of greatest value. 
With instrumentation currently available, an observation 
program concentrated over two months during S2's next 
pericenter passage (2018) will prove to be sufficient as a 
test for the Einstein Equivalence Principle. Figure ([2]) 
argues that a small handful of spectral measurements of 
S2 at 10km s~^ around pericenter imply a la accuracy 
on q; of ~ 0.3. 

The approach we have taken above focuses on the es- 
sentials. The degeneracy between a, I and a for spec- 
troscopy has been lifted by using astrometry only to pro- 
vide /. In practice, all the parameters are fitted simulta- 
neously to both astrometry and spectroscopy. We have 
done simulations to verify that when this is done, the 
degeneracy is implicitly broken by the mechanism high- 
lighted in this Letter. 

Relativistic effects can be expected to become increas- 
ingly important as corrections in other astrophysics re- 
lating to the S stars. Three areas where this can be 
expected are the following. 

1. The combination of spectroscopic and astrometric 
S Star data provides us with the distance to the 
center of the galaxy. Astrometry is sensitive to the 
angular size of the orbit, while spectroscopy on the 
physical size. The quotient is the distance to the 
galactic center (|Eisenhauer et al...2003. ). 



Earth. Spectroscopy therefore has the power to 
determine our velocity with re spect to the central 
black hole (|Angelil et al.|[20Tol) . thus constraining 
the U component of the Galactic local standard of 
rest. 

3. The form of the mass distribution within the inner 
arcsecond affects the S Star orbits. A better under- 
standing of the density profile will provide insight 
into the region's star capture and formation his- 
tory, and to the central dark matter d istribution 
(IGhez et al.l[2008l: iGiUessen et al.ll2009aD . 

In exploring 1 and 2, one cannot easily avoid relativity 
simply by considering stars with larger orbit sizes (3), as 
perturbations due to the enclosed mass become a prob- 
lem. With the accuracy regime that spectroscopy will en- 
ter in the the coming decade, one of two types of pertur- 
bations to the redshift must be faced: either those from 
the extended mass distribution, for S Stars with large or- 
bits, or, those perturbations from relativity, for smaller 
orbit S Stars. While effects due to relativity are well un- 
derstood in principle, and can be easily treated, the con- 
stitution of the extended system is poorly understood, 
and expects a more grueling treatment. In the coming 
decade, spectroscopic S Star accuracy at ~ 10 km s~^ 
is expected to be available, and the discovery of stars 
closer in to the black hole is anticipated. In light of these 
prospects, so that constraints on these quantities may be 
improved, relativistic perturbations, while interesting in 
their own right, can no longer be ignored and must be 
faced. 



2. The position of the observed line depends on the 
velocity of the BH-star system with respect to the 
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